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INTRODUCTION 


The  need  to  reduce  weight  of  structures  has  led  to  the  development  of 
numerous  fiber  reinforced  composite  materials.  Due  to  the  lack  of  a  well 
established  database,  these  new  materials  can  not  be  utilized  to  their  full 
potential.  The  majority  of  available  mechanical  property  data  in  the  open 
literature  was  generated  under  uniaxial  loading  which  represents  a  small 
fraction  of  existing  structural  loading  schemes.  Attempts  to  generate 
composite  property  data  under  biaxial  loading  have  been  made  by  several 
investigators  using  various  specimen  types.  A  review  of  the  multiaxial 
loading  tests  for  composite  materials  is  provided  by  Chen  and  Matthews [1]  .  A 
brief  criticism  of  the  drawbacks  of  these  specimens  was  presented  in  an 
earlier  article^  .  A  new  specimen  for  the  determination  of  the  mechanical 
properties  of  composites  under  biaxial  flexure  loading  was  introduced  in  the 
latter  publications.  Disk  specimens  were  used  to  generate  biaxial  property 
data  for  both  cross-ply [2]  and  unidirectional [3]  aluminum  matrix  composites  and 
a  graphite  fiber/epoxy  cross-ply  composite^.  The  main  advantages  of  the 
disk  specimen  were:  (1)  minimal  machining,  (2)  small  volume  (5  cm  diameter  for 
thicknesses  up  to  3  mm),  (3)  no  special  grips  or  loading  fixtures  required  and 
(4)  near  constant  strain  contours  within  the  loading  area. 

The  disk  specimen  was  used  by  several  investigators  to  conduct  biaxial 
flexure  tests  of  isotropic  ceramic  materials.  These  investigations  included 
analyses  of  the  disk  specimen  by  Finite  Element  Methods  (FEM) .  Three 
variations  of  the  biaxial  flexure  test  exist  for  homogeneous  isotropic 
materials  : 
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Ball-on-Ring  Test 


a . 

b.  Piston-on-Three-Ball  Test 

c.  Ring-on-Ring  Test 

The  piston-on-three-ball  test  is  the  ASTM  ANSI/F394-78  standard  (1978,  part 
43)  .  Recently,  Shetty  et  alt6]  and  Chao  et  al[7]  introduced  the  hydraulic- 
pressure-loading  test  which  loads  the  disk  uniformly. 

The  present  work  evaluates  the  biaxial  strain  field  of  the  disk  for 
continuous  fiber  reinforced  composite  materials  with  both  finite  element  and 
closed  form  analytical  methods.  Effects  of  the  disk  geometry  and  the  material 
characteristics  on  the  resulting  strain  fields  were  explored,  modeled  and 
discussed.  Another  intent  of  this  work  was  to  further  promote  the  adoption  of 
the  disk  and  related  test  procedures  as  a  standard  biaxial  flexure  specimen. 


FINITE  ELEMENT  MODELING 


The  material  considered  for  the  finite  element  analysis  was  a  16-ply 
graphite/epoxy  composite  with  a  balanced  cross-ply  sequence  of  [ ( 02/902 ) 2] s . 

Material  properties  corresponded  to  Hercules  (TM)  3501-6  epoxy  resin/AS-4 
continuous  graphite  fiber  prepreg.  Resin  and  lamina  properties  are  provided 
in  Tables  (1)  and  (2),  respectively.  The  processed  laminate  had  a  60%  fiber 
volume.  Diameter  and  thickness  values  of  the  disk  were  50.8mm  and  2.25mm, 
respectively. 

The  loading  arrangement  for  this  analysis  is  shown  in  Fig.  (1).  Shetty  et 
al[5]  found  that  when  ceramic  disks  deflect  biaxially  under  pin  loading,  the 
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concentric  loading  converted  from  initially  uniform  to  ring-type. 


Ritter  et 


al[8]  made  the  same  conclusion.  Ring-type  loading  was  confirmed  by  the 
appearance  of  ring-shaped  markings  under  the  loading  pin  on  the  compressive 
surfaces  of  tested  diskst2_4].  Thus,  all  strain  profiles  discussed  herein 
corresponded  to  concentric  ring-type  loading. 

Modeling  of  the  16-ply  laminated  disk  was  accomplished  using  the  NISA- 
II[95(TM)  finite  element  analysis  program.  The  Reissner-Mindlin  approach[10] 
revealed  the  problem  of  shear  locking  which  is  associated  with  bending  of  thin 
plates.  Papadopoulos  and  Taylor[11]  described  a  triangular  element  with  a 
unique  interpolational  approach  that  avoids  shear  locking.  The  shear  locking 
problem  was  overcome  in  the  current  model  through  the  use  of  higher-ordered 
element  formulations.  A  combination  of  NISA  type  32,  3-D  laminated  composite, 
quadrilateral  and  triangular  general  shell  elements  were  used.  The 
formulations  used  quadratic  (6-noded  triangles  and  8-noded  quadrilaterals) 
rather  than  linear  interpolation  functions  and  included  transverse  shear 
deformation  capability.  A  ring  of  triangular  elements  was  necessary  to  model 
the  center  of  the  disk  due  to  the  converging  shape  of  the  quadrilateral 
elements  in  this  vicinity.  The  analysis  performed  was  linearly  elastic  and 
incorporated  several  assumptions  from  Classical  Laminated  Plate  Theory^121 
(CLPT) ,  namely:  (1)  the  laminate  consisted  of  multiple  layers  of  perfectly 
bonded  orthotropic  materials,  (2)  no  through-thickness  normal  deformations 
were  allowed,  and  (3)  the  bond  layers  between  plies  were  infinitesimally  thin. 
A  further  assumption  resulting  from  the  NISA- I I  element  formulation  was  that 
shell  elements  possessed  infinite  rotational  stiffness  about  the  normal 
direction  to  their  surface. 


3 


The  disk  specimen  possessed  quarter  symmetry  in  terms  of  material 
properties,  boundary,  and  loading  conditions  collectively.  Therefore,  one 
quadrant  was  discretized  using  330  elements  and  1057  nodal  points.  Global 
axes  (X,Y)  were  taken  parallel  to  the  0°  and  90°  ply  directions,  respectively. 
The  principal  (or  material)  axes  (1,2)  for  each  ply  were  aligned  so  that  the 
1-axis  was  parallel  to  the  fiber  direction  and  the  2-axis  was  perpendicular  to 
the  fiber  direction.  The  displacement  boundary  conditions  are  shown  in  Fig. 
(2) .  In  this  figure,  U,  V,  and  W  represent  the  displacements  along  the  X,  Y, 
and  Z  axes,  respectively  and  U1,  V*,  and  W’  represent  rotations  along  the  same 
axes .  Quarter  symmetry  was  enforced  by  restraining  the  following  Degrees  of 
Freedom  (DOF) : 

DOF  LOCATION 
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at 
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II 
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II 
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II 
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II 
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o 

II 
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The  outer  edge  of  the  disk  was  restrained  from  displacement  in  the  Z- 
direction  by  setting  W=0  at  nodes  along  this  edge.  A  distributed  load  of 
4,45kN  was  applied  as  a  series  of  concentrated  forces  at  nodal  points  defining 
a  circle  of  radius  "b".  Radius  "b"  was  equal  to  the  load  pin  radius  (6.1mm) 
used  in  tests  by  Tsangarakis  et  al[2~4J.  In  further  accordance  with  CLPT[12], 
the  in-plane  strains  were  assumed  linearly  dependent  with  the  distance  from 
the  neutral  plane.  The  neutral  plane  was  assumed  to  be  coincident  with  the 
mid-plane  of  the  disk  according  to  CLPT  for  symmetric  laminates. 

The  deflected  disk  shape,  superimposed  with  the  unloaded  configuration,  and 
the  transverse  (out-of-plane)  displacement  contours  are  shown  in  Figs.  (3)  and 
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(4),  respectively.  Note  the  asymmetrical  profile  of  the  transverse 
displacement  contours  in  Fig.  (4)  which  resulted  directly  from  the  anisotropy 
of  the  laminate.  Figure  (5)  shows  in-plane  stress  contours  for  the 
compressive  surface  layer.  The  compressive  surface  layer  alone  was  examined 
since  stress  and  strain  fields  of  the  remaining  layers  varied  linearly  with 
the  distance  from  the  laminate  neutral  plane.  Stress  output  was  automatically 
extrapolated  from  the  Gauss  (integration)  points  to  the  nodes,  while  strains 
were  obtained  at  element  centroids.  In-plane  strain  versus  radial  distance 
plots  of  Figs.  (6-8)  show  principal  and  shear  strains  plotted  against  rays  of 

element  centroidal  locations.  Ray  angles  were  measured  from  the  global  X- 
axis.  Results  indicated  that  principal  strains  e11  =  -0.0137  m/m  and  e22  = 

-0.0155  m/m  and  formed  a  nearly  equi-biaxial  strain  field  within  the  load  pin 
radius.  At  the  center  of  the  disk,  the  two  principal  strains  varied  by  12%. 
The  in-plane  shear  strain  at  any  orientation  within  the  load  pin  radius 
converged  to  zero.  Localized  effects  of  the  strain  field  due  to  the  applied 
loading  appeared  as  knees  on  each  strain  component  curve  and  decayed  towards 
the  center  approaching  constant  values.  Furthermore,  the  principal  strains 
became  invariant  to  angular  position  for  points  approaching  the  center. 
Finite  element  modeling  by  Ritter  et  alt8J  and  experimental  tests  by  Shetty  et 
al^l  confirmed  the  existence  of  similar  localized  stress  magnifications  at 
the  load  point  for  isotropic  ceramic  disks.  Outside  of  the  load  pin  radius, 
the  principal  and  shear  strains  had  significant  dependence  upon  angular 
position.  These  strain  profiles  were  in  agreement  with  strain  gage 
measurements  made  Tsangarakis  et  al^2_4].  The  peak  out-of-plane  deflection  for 
a  load  of  4.45kN  was  2.835xl0~3  m  occurring  at  the  center  of  the  disk. 
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A  special  case  of  interest  was  investigated  in  which  a  symmetric  cross-ply 
laminate  stacking  sequence  was  optimized  to  provide  an  equi-biaxial  strain 
field  within  the  load  pin  radius.  In  particular,  the  stacking  sequence 
resulting  in  equivalent  laminate  flexure  moduli  ELF  and  ETF  (where  ELF 
designated  the  longitudinal  or  X-direction  flexure  modulus  and  ETF  designated 

the  transverse  or  Y-direction  flexure  modulus)  was  determined.  For  symmetric 
laminates,  the  following  compliance  equations  by  Tsai^13^  were  used  to 
determine  the  flexure  moduli  ELF  and  ETF : 

Elf  =  l/dn*  and  ETF  =  l/d99* 

[d*]  =  [d]h3/12  =  Normalized  Bending  Compliance 
h  =  Laminate  Thickness 

[d]  =  [D]”1  =  Laminate  Bending  Compliance  Matrix 
[D]  =  Laminate  Bending  Stiffness  Matrix 

Figure  (9)  shows  the  variations  in  flexure  moduli  ELF  and  ETF  for  4  different 
symmetric  cross-ply  laminates.  The  [  (0/90) s/ ( 90/0) s]  s  laminate  yielded 
equivalent  values  of  ELF  and  ETF  at  63.30  GPa .  The  previous  FEA  model  was 

updated  to  include  the  optimized  ply  stacking  sequence  and  re-executed. 
Compressive  surface  principal  strains  sn  and  e2?  were  equal  to  -0.0143  m/m 

within  the  load  pin  radius  and  became  invariant  to  angular  position  for  points 
approaching  the  disk  center.  Peak  out-of-plane  deflection  was  2.805xl0~3  m 
occurring  at  the  center  of  the  disk. 


where : 


and: 
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MODELING  THE  EFFECT  OF  FIBER  SIZE  AND  SPACING 


Considering  a  0°  ply,  the  two  neighboring  fibers  shown  in  Fig.  (10)  have 
radius  R  and  are  separated  by  matrix  material  of  thickness  h.  The  coordinate 

system  of  this  figure  is  the  same  as  that  used  in  the  finite  element  model 

(see  Fig.  2)  .  Let  the  fibers  be  designated  as  i  and  i+1  and  experience  the 
strains  e^tx^zj  and  sxx  (x,  y,  zi+1)  ,  respectively.  Further  assume  that  the: 

1.  Fibers  are  well  bonded  to  the  matrix. 

2.  Fibers  and  matrix  obey  Hooke* s  Law. 

3.  Deformation  "w"  along  the  z-axis  is  negligible  compared 

to  deformations  "u"  and  "v"  along  axes  x  and  y,  respectively. 

4.  Principal  strains  on  the  tensile  side  of  the  disk  are  equal 

and  constant  on  a  given  plane  and  within  the  load  pin  region. 

5.  Fibers  are  subjected  to  only  uniform  axial  loads. 

The  interlaminar  shear  strain  yzx  experienced  by  the  matrix  will  be: 

yzx  =  dw/dx  +  du/dz  (1) 

From  the  third  assumption,  the  displacement  w  and  its  derivative  along  the  z- 
axis  are  negligible.  Thus,  the  first  term  in  the  right  side  of  Eq.(l)  may  be 
omitted.  Using  the  model  discussed  by  Puppo  and  Evensen[14],  Garg  et  al[15], 
and  recently  by  Gao  and  Reif snider [16]  for  a  balanced  laminated  composite,  the 
interlaminar  shear  strain  can  be  approximated  as: 

Yzx  =  du/dz  =  (ui+1-ui)/h  (2) 


7 


Let  zi  be  the  distance  of  fiber  i  from  the  mid-plane. 


If  t  is  the  disk 


thickness  and  sXX(t/2)  is  principal  strain  of  the  outermost  fiber  ( 

then : 

sxx(x,y,Zi)  =  (2zi/t)  exx(t/2) 

The  strain  e  (x,y,  may  also  be  expressed  as: 

exx  (x,  y,  Zi)  =  dui/ dx 

Integrating  Eq. (4) : 

=J  8;<x(x,y,zi)dx 

O 

Substituting  the  principal  strain  from  Eq. (3)  into  Eq. (5) : 

ui  =J  sxx(t/2)  (2zx/t)dx 

J  O 

And  because  sxx(t/2)  is  constant  per  assumption  4: 

Ui  =  Sxx(t/2)  I2zt/t)x 
Ui  +  1  =  exx(t/2)  (2zi  +  1/t)K 

Substituting  Eqs.(7a,b)  into  Eq.{2): 

Yzx  =  2hexx(t/2)  Izi+l-zi]x/t: 

From  Fig. (10),  the  difference  (zi  +  1  -  zt)  may  be  expressed  as  (2R  +  h) . 
Substituting  into  Eq. (8) : 

Yzx  =  2exx(t/2)  [i+(2R/h)]x/t 


z=t/2 )  , 

(3) 

(4) 

(5) 

(6) 

(7a) 

(7b) 

(8) 
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The  same  consideration  is  followed  for  shear  strain  yzy  in  a  90°  layer  where 
the  x  subscripts  and  coordinate  term  are  replaced  by  y.  This  treatment,  of 
course,  is  valid  for  balanced  composites. 

The  shear  strains  yzx  and  Yzy  in  90°  and  °°  laYers'  respectively,  are  found 
through  the  usual  treatment  with  the  theory  of  bending  of  thin  composite 
plates [15].  These  strains  are  maximized  at  locations  designated  as  points  A 
in  Fig.  (11),  that  is,  where  the  matrix  thickness  h  is  minimal.  At  remote 
locations  designated  as  points  B,  the  shear  strain  is  minimal. 

Equation  (9)  indicates  that  the  interlaminar  shear  strain  depends  on  the 
disk  thickness,  fiber  radius,  and  interfiber  matrix  thickness.  The 
interlaminar  shear  strain  also  depends  on  distance  x  and  the  strain  eXX(t/2)* 

For  failure  to  occur  by  matrix  shearing,  the  interlaminar  shear  strain  must 
exceed  a  critical  value  y* .  This  is  the  value  of  the  shear  strain  required  to 
form  a  delamination.  From  Eq.(9),  it  is  evident  that  as  the  ratio  2R/h 
increases,  the  matrix  shear  strain  increases.  Note  that  this  shear  strain  is 
independent  of  the  distance  from  the  fibers  to  the  neutral  plane.  Larger 
fibers  will  favor  the  formation  of  delaminations.  When  fibers  are  bent  or 
misaligned,  the  matrix  layer  thickness  h  may  decrease  between  two  neighboring 
fibers  leading  to  delamination.  Thus  the  ratio  2R/h  may  be  considered  as  a 
local  shear  strain  concentration  factor,  KLS.  Because  this  ratio  may  vary  at 

many  locations  within  the  disk,  formation  of  multiple  delaminations  is 
possible. 
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DISCUSSION 


The  strain  profiles  depicted  in  Figs. (6,7)  were  verified  by  Finite  Element 
Analysis  and  experimentally  with  strain  gages [2' 4 ]  for  a  balanced,  cross-ply 
[(02/902)2]s  fiber  reinforced  composite.  The  two  principal  strains  e11  and  g22 

were  found  to  vary  within  12%  and  were  constant  within  the  load  pin  region. 
These  observations  were  also  found  for  isotropic  ceramic  disks  by  Ritter  et 
al  ^  .  The  two  principal  strains  were  zero  at  the  disk  mid-thickness  .  The 
ability  to  tailor  fiber  reinforced  composite  laminates  resulted  in  an 
optimized  ply  stacking  sequence,  namely  [( 0/90 ) s/ ( 90/0) s]  s,  which  yielded  an 

equi-biaxial  principal  strain  field  within  the  load  pin  radius. 

In  the  unidirectional  aluminum  composite  disk[3^,  the  maximum  principal 
tensile  strain  (parallel  to  the  fibers)  was  constant  within  the  load  pin  area 
as  observed  by  strain  gage  measurements.  The  principal  tensile  strain 
perpendicular  to  the  fibers  was  maximal  at  the  disk  center  and  decayed  to  zero 
thereafter  towards  the  disk  edge. 

As  Eq.(9)  indicates,  within  the  load  pin  region,  the  interlaminar  shear 
strain  increases  with  distance  from  the  center  of  the  disk.  Letting  x  =  b, 
where  b  is  the  load  pin  radius,  the  interlaminar  shear  strain  at  the  edge  of 
the  load  pin  is: 


Vzx  =  sxx(t/2)  d  +  2R/h)2b/t 


(10) 


The  interlaminar  shear  strain  is  maximized  at  the  disk  edge  and  is 
approximated  as: 
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(1  +  2R/h) ( r-b) / t 


(11) 


Yzx 


2e 


xx ( t/2 ) 


Since  the  interlaminar  shear  strain  is  dependent  upon  the  ratio  2R/h,  large 
or  closely  spaced  fibers  will  induce  higher  shear  strain  values •  Similar 
conclusions  were  drawn  by  Ramsteiner  and  Theysohn[17]  for  a  short  glass  fiber 
reinforced  organic  composite.  These  investigators  found  that  as  the  fiber 
radius  increased,  the  failure  tensile  strength  decreased.  Thus,  evidence 
exists  that  larger  fibers  entail  higher  local  shear  strain  concentrations 
(Ktc).  The  ratio  2R/h  depends  on  the  fiber  volume  fraction,  Vf,  of  the 

composite.  A  simple  geometrical  consideration  shows  that  the  ratio  dependence 
on  V  is: 


2R/h  =  2/(/tc7v7  "2)  (12) 

The  above  equation  indicates  that  as  the  fiber  volume  fraction  increases, 
so  does  the  aspect  ratio  2R/h.  For  fiber  volume  fractions  of  0.1  and  0.6,  the 
respective  values  of  the  aspect  ratio  are  0.555  and  6.940.  These  values 
represent  the  lower  and  upper  bounds,  respectively,  of  the  local  shear  strain 
concentration  factor  KLS  which  is  plotted  in  Fig.  (12)  . 

It  was  shown  that  the  interlaminar  shear  strain  is  independent  of  the 
distance  from  the  disk  mid-plane.  When  the  value  of  the  interlaminar  shear 
strain,  as  given  by  Eq.(9),  exceeds  a  critical  value,  a  matrix  crack  may  form 
by  shear.  Because  KLS  may  vary  greatly  (0.55  to  6.940)  due  to  uneven 

distribution  of  the  fibers,  the  matrix  shear  strength  may  be  exceeded  at  many 
locations  simultaneously.  Thus,  multiple  cracks  may  form  at  different 
locations.  Such  cracks  are  shown  in  Fig.  (13)  for  a  silicon  carbide  fiber 
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(0.14mm  diameter)  reinforced  aluminum  composite^,  and  in  Fig.  (14)  for  a 
carbon  fiber  reinforced  resin  composite^.  In  Fig.  (15)  ,  a  C-scan  of  a  failed 
organic  composite  disk  shows  the  multiplicity  of  the  generated  matrix  cracks. 
Existing  theories  as  those  of  mixed  modes  fracture  energy,  shear  lag,  and  the 
empirical  mixed  mode  fracture  toughness  of  Wang  et  al^18^  do  not  explain  the 
formation  of  multiple  delaminations. 

An  estimate  of  the  interlaminar  shear  strength  of  the  organic  composite^ 
was  obtained  by  testing  two  sets  of  short  beam  shear  specimens  in  accordance 
with  the  ASTM  D2344  procedure.  The  specimen  size  was  1.9mm  x  6.3mm  x  25mm. 
Each  set  consisted  of  six  coupons.  The  estimates  of  the  interlaminar  shear 
strength  were  44.5  +  9.5  MPa  and  38.7  +  7.5  MPa  for  the  first  and  second  set, 
respectively.  The  average  was  41.6  MPa  which  represented  a  gross  estimate  of 
the  interlaminar  shear  strength  of  the  composite. 


CONCLUSIONS  AND  RECOMMENDATIONS 


The  strain  fields  in  fiber  reinforced  laminated  composite  disks  loaded  with 
a  pin  on  one  side,  while  the  other  side  was  simply  supported  cicumf erentially, 
were  examined.  Finite  element  analysis  results  indicated  that  in  a  balanced, 
[(02/902)2]s  cross-ply  composite,  the  in-plane  principal  strains  were  nearly 

equal  and  constant  within  the  load  pin  radius.  An  optimized  cross-ply 
stacking  sequence  of  [  (0/90) s/ (90/0) s]  s  resulted  in  an  equi-biaxial  principal 

strain  field  within  the  load  pin  radius.  The  interlaminar  shear  strains  were 
shown  to  be  constant  through  the  disk  thickness  and  greatly  dependent  upon 
interfiber  distance  and  fiber  diameter.  It  was  shown  that  large,  closely 
spaced  fibers  could  intensify  the  shear  strain  field  to  the  extent  of  forming 


12 


single  or  multiple  delaminations.  Because  perfect  fiber  alignment  can  not  be 
achieved  throughout  the  laminate  thickness,  the  present  test  method  can  not  be 
used  to  determine  the  interlaminar  shear  strength  of  fiber  reinforced 
composites .  However,  the  disk  specimen  would  be  useful  in  comparative  studies 
as  a  quality  control  tool  and  to  assess  the  response  of  the  composite  to 
biaxial  flexure  loading. 
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TABLE  1 


PROPERTIES  OF  THE  3501-6  RESIN 


TABLE  2 


Ultimate  Tensile  Strength,  UTS  =  70  MPa 
Elastic  Modulus,  E  =  4.4  GPa 
Tensile  Strain-to-Failure,  sf  =  0.017  m/m 


-  UNIDIRECTIONAL  PROPERTIES  OF  3501-6  EP/60%  AS-4  CARBON  FIBER  PREPREG 

Elastic  Modulus,  E11  =  115  GPa 
Elastic  Modulus,  E22  =  11  GPa  [15^ 


Major 

Poisson  1 s 

Ratios , 

v12  _  V1 3  ~  0.27 

[16] 

Minor 

Poisson ' s 

Ratios , 

V21  =  V31  =  0.02 

[16] 

Ultimate  Longitudinal  Tensile  Strain,  su  =  0.0156  m/m  [161 
Ultimate  In-Plane  Shear  Strain,  s12  =  0.0057  m/m  [161 
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Fig. (3)  DEFORMED  AND  UNLOADED  DISK  GEOMETRIES 
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Fig. (4)  TRANSVERSE  DISPLACEMENTS  OF  DISK  MID-PLANE 
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Fig.  (9)  FLEXURE  MODULI  vs.  CROSS-PLY  LAMINA  STACKING  SEQUENCE 


(Z-AXIS  PARALLEL  TO  DISK  THICKNESS) 


Fig.  (10)  EFFECT  OF  FIBER  DIAMETER  AND  SPACING 
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Fig.  (11)  LOCATIONS  OF  MAXIMUM/MINIMUM  LOCAL  SHEAR  STRAINS 


LOCAL  SHEAR  STRAIN  CONCENTRATION  K  vs.  FIBER  VOLUME  V 


Fig.  (13)  INTERLAMINAR  SHEAR  CRACKS  FORMED  IN  SILICON 
CARBIDE /ALUMINUM  COMPOSITE  DISK 


Fig.  (14)  INTERLAMINAR  SHEAR  CRACK  FORMED  IN 
GRAPHITE /EPOXY  COMPOSITE  DISK 


Fig. (15)  C-SCAN  OF  GRAPHITE /EPOXY  DISK  WITH  MULTIPLE  CRACKS 
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